Interaction-enhanced ferromagnetic insulating state of the edge 
of a two-dimensional topological insulator 
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The spin-related aspects of electron-electron interactions at the edge of a two-dimensional topo- 
logical insulator are addressed. In the absence of the magnetic field, the interacting edge electrons 
form a helical Luttinger liquid with gap less collective excitations that carry both spin and charge. 
The edge exhibits in-plane XY ferromagnetic correlations but the long-range order is not formed 
and the edge is a perfect conductor. Applying the magnetic field orders the edge ferromagnetically 
and opens a gap A in the spectrum of collective spin-charge excitations, which can be substantially 
enhanced by the interactions. The exponent of the scaling dependence A oc H 1 ^ 2 ~ K ^ > of the gap on 
the magnetic field H is controlled by the Luttinger liquid interaction parameter K, which can be 
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extracted from the activation behavior G xx ~ ^^exp(— A/T) of the longitudinal conductance at 
lower temperatures T< A. This provides a simple way to probe the interactions experimentally. 

PACS numbers: 72.25.-b, 71.10.Pm, 73.43.Lp 
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Introduction. Topological insulators fll-[l5| form a new 
class of materials with nontrivial band structure caused 
by spin-orbit interactions. The key feature that distin- 
guishes a topological insulator from a conventional one 
is the presence of gapless single-particle surface or edge 
states. The edge of a two-dimensional (2D) topologi- 
cal insulator [Jo, 0? Is| supports two branches of gap- 
less counter-propagating helical states with opposite spin 
projections on the axis perpendicular to the plane of the 
sample (Fig. [I]). Protected by the time-reversal symme- 
try (TRS) against single-particle backscattering [la, Il7j. 
these edge modes serve as nearly ideal conducting chan- 
nels that give rise to the quantum spin Hall (QSH) ef- 
fect. So far, the 2D topological insulator was realized 
in HgTe-CdTe quantum wells, which was first predicted 
theoretically |9[ and shortly after confirmed experimen- 
tally 0, Ell • 

The interactions between electrons in the counter- 
propagating states lead to a novel helical Luttinger liquid 
(LL) phase 
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23j with unusual spin-related behavior. 
It would be beneficial to detect the signatures of interac- 
tions in the experiment. However, in the critical LL state 
the collective excitations are gapless and the edge con- 
ductance remains essentially unaffected by the interac- 
tions [24j,[25|. Therefore, as long as TRS is preserved, the 
interactions do not reveal themselves in the longitudinal 
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conductance G xx = 2 j-^ and probing them by transport 
measurement requires creating a tunneling setup, such as 
a point contact between two QSH edges [20l422j. 

The situation changes if one breaks TRS. Experimen- 
tally, this may easily be achieved by applying the external 
magnetic field [To|, [H|. The magnetic field couples the 
helical states, spin polarizes the edge in the plane of the 
sample, and opens a gap in the single-particle spectrum. 
As will be shown here, in the absence of the magnetic field 
the interactions in the helical liquid result in the tendency 
towards in-plane XY ferromagnetism manifested in the 




FIG. 1: Helical edge states of a 2D topological insulator. 
(Left) The states propagating in the opposite directions have 
opposite spin projections on the direction perpendicular to 
the plane of the sample . (Right) In the absence of the mag- 
netic field the counter-propagating states are gapless. Shaded 
regions depict the continuum of the extended bulk states with 
the insulator gap eo- 



power-law decay of the spin correlations. Therefore, elec- 
tron interactions will naturally facilitate the formation 
of the ferromagnetic state once the magnetic field is ap- 
plied. This is revealed in the enhancement of the gap in 
the spectrum of collective excitations as compared to the 
noninter acting case. The gap leads to the insulating be- 
havior of the edge at low temperatures. The suppression 
of the longitudinal conductance with the applied mag- 
netic field was indeed observed experimentally in HgTe 
quantum wells Io|, 11]. 

The present paper is devoted to the theory of this ferro- 
magnetic insulating state of the edge of a 2D topological 
insulator, induced by the magnetic field and enhanced 
by the interactions. The obtained scaling dependence 



of the transport gap A ex H 1 ^ 2- ^ on the magnetic 
field H allows one to directly extract the Luttinger liq- 
uid interaction parameter K from the activation behavior 
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G xx ~ 7^exp(— A/T) of the longitudinal conductance. 
A device in the Hall-bar geometry with the applied mag- 
netic field thus provides a minimal setup to probe the in- 
teractions and the predicted dependence could be tested 
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already on the existing data ToL Hl|. 

Model and Hamiltonian. The effective low-energy 
Hamiltonian for the interacting edge electrons in the 
presence of the magnetic field [l8| may be written 
down in the helical basis of right-moving (with respect 
to the x direction along the edge) spin-up (t) and left- 
moving spin-down (\) states as 



H = H() + Hrr 



H m = -A J 



H[, Hq 



dxip^(x)vpa z i/j(x), (1) 



dxip^ (x)(a x coscpo + cF y sin (fo)ip(x), (2) 



Hi = \J dxdx'^l(x)^l{x')V{x 



■ x')^ af (x')^ a (x). 

(3) 

Here, -0* = (ip^^ipi) is the two-component fermionic field 
operator, p = —ihd x , and o- Xl a yi a z are the Pauli ma- 
trices in the helical basis. The part H m describes the 
effect of the external magnetic field. For the in-plane 
orientation, H = H(cos cpo, sin y?o? 0)? on ly the Zeeman 
effect is present, whereas the orbital effect vanishes; the 
angle (po correspond to the direction of the field in the 
xy plane of the 2D sample and the gap is given by the 
Zeeman energy A n ~ HbH. In case of the perpendic- 
ular orientation of the field, H = (0,0,i^), the Zeeman 
effect does not affect the dynamics and only the orbital 
effect remains. The orbital effect of the perpendicular 
field is estimated [11] to be stronger than the in-plane 
Zeeman effect, Ao± ~ lOAoy; A n ~ 3K and Ao± « 30K 
at i7 = IT. For arbitrary field-orientation, the single- 
particle gap Ao scales linearly with the magnetic field, 
A oc H. 

We consider the case of Coulomb interactions, V(x) = 
el/\x\ in Eq. (J3j), possibly screened by the nearby metallic 



electrodes beyond some length l s \ the charge e* = e/y/H 
incorporates the effects of screening by the dielectric envi- 
ronment. This allows us to consider both unscreened and 
screened interactions, the latter modeling practically any 
finite-range interactions. The short-scale spatial cutoff a 
of the theory [Eqs.(pQ), ((2]), and ((3])] and of the poten- 
tial V(x) is set by the decay scale of the edge states into 
the bulk. For simplicity, it is assumed that the chemical 
potential is exactly at the branch crossing e = of the 
unperturbed edge spectrum e p = ±vp, where the correla- 
tion effects are strongest. This can be achieved by tuning 
the gate voltage to the minimum of the longitudinal con- 
ductance. 

The Hamiltonian H [Eqs. (PQ), ©, and (j3])] describes 
one-dimensional interacting Dirac fermions, which are 
massive in the presence of the magnetic field. For point 
interactions, this is known as the Thirring model [28|,|29j. 
This model can be mapped to the sine-Gordon model by 
mean of the bosonization procedure [28, 29]. One relates 




FIG. 2: Collective spin-charge excitations of the edge of a 
2D topological insulator. Excitations are described by the 
phase variable <p(x), which determines both the in-plane spin 
polarization [Eq. J5J)] and charge density [Eq. J5J)]. A kink of 
height TT in tp(x) rotates the spin polarization in the xy plane 
of the sample by 2ty and accumulates a unit charge in the 
region of variation of ip(x). 



the fermion fields xp^^x) of the right and left movers to 
the bosonic ones ip^^(x) as 



\j2ixa 



(4) 



where the Klein factors are omitted. The operators 
tp(x) = \[^{x) + <Pi(x)] and 0(x) = \[^{x) - ip±{x)\ 
satisfy the canonical (up to a coefficient) commutation re- 
lations [<p(x), d x '0(x')] = — \tt5(x — x') and are related to 
the coordinate and momentum variables of the collective 
excitations. In terms of (p(x) and 0(x), the Hamiltonian 
H [Eqs. (P), (j2j), and (j3j)] can be expressed as 



h = A / dxv[(d x e) 2 + (d xV ) 2 



Hm = — - 



dx cos[2(p(x) + <£o]> 



(5) 



(6) 



Hi = ^ J dxdx' dMx)V(x - x')d x ^(x'). (7) 

The sine-Gordon Hamiltonian, Eqs. ©, ©, and (|6j), de- 
scribes the dynamics of the collective edge excitations of 
a 2D topological insulator in the presence of the mag- 
netic field. Further, we proceed with the analysis of this 
Hamiltonian. 

Collective spin- charge excitations. To visualize the col- 
lective excitations described by Eqs. (j5j), (j6j), and (0, let 
us link the fields (f(x) and Q(x) to the physical observ- 
ables. From the relation (j4]), one obtains 



Sx(%) \ = 1 / cos(-2ip(x)) 
s y( x ) J ~ Z^a \ sm(-2ip(x)) 



(8) 



for the x and y components of the spin density operator 
s(x) = ipl(x)cr CTCT filj CT f(x) (defined without 1/2 factor) and 

s z (x) = -d x 0(x), p(x) = -d x <p(x) (9) 

7T 7T 

for the z component of the spin density and the parti- 
cle density p(x) = ^J(x)^ cr (x) operators. As seen from 
Eq. (|8|), the angle — 2tp(x) corresponds to the direction of 
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the spin polarization in the xy plane and the field ip(x) is 
thus directly related to the spin degrees of freedom. At 
the same time, according to Eq. (|9]), the charge density is 
determined by the gradient of <p(x). Therefore, the col- 
lective excitations carry both charge and spin, which is a 
direct consequence of the coupling between the spin and 
orbital degrees of freedom in the single-particle states. 
A kink of height ir in cp(x) rotates the spin polarization 
in the xy plane by 2tt and simultaneously accumulates a 
unit charge in the region of variation of <p(x), Fig. [2] It 
was suggested in Ref. [18| to exploit this bonding of spin 
and charge degrees of freedom to observe charge fraction- 
alization effects in domain-wall structures with inhomo- 
geneous magnetization. 

Gapless helical liquid, H = 0. Let us first consider the 
system in the absence of the magnetic field, H m = 0, 
when the edge is in the helical LL phase, and obtain the 
excitation spectrum and basic correlations. The calcula- 
tions can be conveniently performed in the Langrange 
finite-temperature formalism. From Eqs. (0 and (0, 
the action for the Fourier transformation ip(u n ,q) = 

J VT dr / dxe iUnT - iqx (p(r,x) (hw n = 2irTn, n e Z) of 
the phase field takes the form 

S [<p] + S i [<p] = Tj2 



2-KK a 



(10) 

The momentum-dependent velocity u q and LL interac- 
tion parameter K q are given by 



u q /v = 1/K q = y/l + V(q)/(irhv) = y/r s In[l/(g,a,)], 

(11) 

where V(q) — 2e 2 ln[l/(g*a)] is the Fourier transform 
of the potential V(x), r s = 2e\j{i\hv) is the Coulomb 
parameter, q* = max(|g|, 1/Z S ), and a* ~ ae _1 / rs . 

From Eqs. ([TP]) and ([TT]) . one obtains the excitation 
spectrum uj(q) = u q \q\ of the collective edge excitations 
of a 2D topological insulator. For unscreened Coulomb 
interactions V(q) = 2e\ ln[l/(|g|a)] at ql s > 1, u q and 
K q depend logarithmically on q and the excitations have 
a ID plasmon-type spectrum uj(q) oc q^/hi(l/q). At spa- 
tial scales exceeding the screening length l S: ql s < 1, 
the interactions become effectively short-range with V(q) 
saturating to the value V(q < = 2e^ ]n(l s /a). The 

velocity u q = u and interaction parameter K q = K be- 
come ^-independent, u/v = 1/K = ^/r s ln(Z s /a*), and 
the spectrum cj(q) = u\q\ linear. In the absence of the 
magnetic field the spectrum is gapless, but for unscreened 
Coulomb interactions the log-dependence of the velocity 
u q signals of a strong tendency towards gap opening. 

Let us now study the correlations. The operators 
that describe coupling between the counter-propagating 
helical modes are given by the "spin-flip" components 
s±(x) — s x (x) ± is y (x) of the spin density (j8j), 

s+(x) = ^\{x)^(x) = "2^-- (12) 




FIG. 3: Ferromagnetic insulating state of the edge of a 2D 
topological insulator. Magnetic field couples the spin up and 
down helical states, opens a gap Ao in the single-particle edge 
spectrum (right), and orders the edge ferromagnetically in the 
plane of the sample (left). The many-body gap A [Eq. JT4J)] 
is enhanced by the interactions compared to the bare gap Ao. 



The pairing tendency is thus directly related to the spin 
polarization in the xy plane of the 2D sample. Calculat- 
ing the correlation function of s±(x) with respect to the 
action (JTJ5J) at zero temperature T = 0, we obtain 



( S+ (x) S _(0))oc 



exp 



<l. 



-4y/]n(\x\/a*) /r 8 

% (ls/\x\f K , \X\>1 S . 

(13) 

For screened Coulomb interactions at \x\ > l s the cor- 
relations ([T3]) of the in-plane spin density s x ^ y (x) have a 
LL power-law decay. For unscreened Coulomb interac- 
tions at \x\ < l S: the decay is slower than any power law. 
The interactions in the helical liquid thus result in the 
tendency towards in-plane XY ferromagnetic ordering. 
However, due to strong quantum fluctuations in a ID sys- 
tem the long-range order is not formed, (s(x)} = 0. For 
unscreened Coulomb interactions, the tendency towards 
ferromagnetism is as strong as that towards Wigner crys- 
tallization in a conventional one-dimensional electron 
system [26|, [27J . Note that numerical factors in the spec- 



trum uj(q) 



•i q \q\ [Eq. ([TT]) ] and correlation function 
([T3]) differ from those of Refs. {if! because in our 
case electrons are single-flavored. 

In the massless LL phase, the edge conductance G = 

2 

is essentially unaffected by the interactions and 
the edge remains a perfect conducting channel [24], EBJ . 
Therefore, in the absence of perturbations that break 
TRS, the interactions do not reveal themselves in the 
transport measurement of the longitudinal conductance 
G xx = 2G of the sample (factor 2 due to two edges). 

Gapped ferromagnetic phase, H > 0. The situation 
changes, if the magnetic field is applied, Ao > in 
H m [Eq. fl5)]. Even in the absence of interactions, the 
magnetic field couples the helical counter-propagating 
states [ll|, LL8| according to Eq. ([2]) and opens a gap Ao in 
the single-particle spectrum e p = d=y (vp) 2 + Aq of the 
Hamiltonian i^o + H m , Fig. [3] In the ground state, the 
edge becomes spin polarized in the plane of the sample 
in the direction </?o, ( s (^)) oc (cos <po , sin (fo , 0) . 
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Opening of the single-particle gap A has a direct 
consequence on transport. For the noninter acting elec- 
trons, the edge conductance can be calculated using the 
Landauer formula and for long enough edge of length 
L ^> hv/ Aq it is given by 



v J 27rfiexp(A /T) + l' 

The presence of the gap makes the edge insulating at tem- 
peratures T « Ao, where the conductance follows the 
Arrhenius activation law G(T) « (^§^j 2exp(— A /T). 

Let us now take the interactions into account. In terms 
of the collective excitations, the effect of the magnetic 
field is described by the cosine term ((6]) in the bosonized 
Hamiltonian. The fact that the ground state is spin po- 
larized means that the phase field ip(x) is locked in the 
minimum of the cosine term, ((p(x)) = —cpo/2. The col- 
lective excitations are now massive and for low energies 
described by the fluctuations of ip(x) around this mini- 
mum. Since even without the magnetic field the interac- 
tions tend to order the edge ferromagnetically, naturally, 
the gap A in the spectrum of the collective excitations 
turns out to be enhanced compared to its bare single- 
particle value Ao. For screened Coulomb interactions we 
obtain 

/ An\^ W 1 

A-eo(-^J cx H— . (14) 

up to a numerical factor ~ 1. Here eo is the bulk insulator 
gap, which determines the high energy cutoff of the edge 
spectrum and is assumed eo ^> Aq. For HgTe quantum 
wells, it is estimated e ~ 100K [11]. The result ([T4j) 
can be obtained by several means, e^., using the self- 
consistent harmonic approximation [29|. 

The gap ([T4]) has a power-law dependence on the bare 
gap Ao ~ HbH and hence on the magnetic field H. The 
exponent 1/(2 — K) of this dependence is controlled by 
the LL interaction parameter K, which varies between 
K = 1 in the noninteracting case and K = for infinitely 
strong finite-range interactions; these cases give the low- 
est A min = A and highest A max ~ y/ Aoeo y/H pos- 
sible values of the many-body gap A, respectively. Due 
to the long-range nature of the Coulomb forces, for un- 
screened interactions the gap appears to be close to A max 
even for moderate interaction strength r s ~ 1. Perform- 
ing the harmonic approximation [29|, we obtain 

A 2 ~ A e exp[-V21n(eo/A )/r a ]. (15) 

The gap (|15p differs from the K = limit A max of 
Eq. ([14]) only by a function of Ao/eo that varies slower 
than any power law. The result (fT5j) applies if the corre- 
lation length /a = hv/A determined from Eq. (fT5j) does 
not exceed the screening length, I a < l s . Otherwise, 
what concerns the gap, the interactions are effectively 



screened and the gap is given by Eq. (fT4|) . For unscreened 
Coulomb interactions, the enhancement of the gap could 
thus be quite substantial: for Ao ~ IK and eo ~ 100K 
one gets A 

max 

10K. The enhancement of the gap 
means, in particular, that interactions should favor ob- 
servation of the effects predicted in Ref. [18|. 

Experimental manifestation. The main practical con- 
sequence of the interactions is that the many-body gap A 
[Eq. JT2J)] plays the role of the transport gap and enters 
the activation law G XX (T) oc 7^ exp(— A/T) for the lon- 
gitudinal conductance at temperatures T <C A. There- 
fore, the scaling dependence (fT4|) of the gap A on the 
magnetic H should be directly accessible in the transport 
measurement allowing one to extract the Luttinger liquid 
interaction parameter K. This provides a direct way to 
probe the interactions in the transport experiment and 
the predicted dependence could be tested already on the 
existing data [Io|, [H| • 

Author is thankful to Konstantin Matveev for valuable 
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under Contract No. DE-AC02-06CH11357. 
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